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ABSTRACT 

The iron Ka line is perhaps the most important spectral diagnostic available 
in the study of accreting black holes. The line is thought to result from the 
reprocessing of external X-rays by the surface of the accretion disk. However, 
as is observed in the solar corona, illumination by energetic particles may also 
produce line emission. In principle, such a process may be uncorrelated with the 
observed X-rays and could explain some of the unexpected variability behavior 
of the Fe Ka line. This paper compares predictions of iron Ka flux generated 
by impacting electrons and protons to that from photoionization. Non-thermal 
power-laws of electrons are considered as well as thermal distributions of electrons 
and virialized protons. The electrons are thought to originate in a magnetically 
dominated accretion disk corona, while the protons are considered in the context 
of a two phase (hot/cold) accretion scenario. In each case, the Fe Ka flux from 
particle impact is found to be < 1% of that produced by photoionization by 
a hard X-ray power-law (normalized to the same energy flux as the particles). 
Thus, the electrons or protons must strike the disk with 10^-10^ times more 
energy flux than radiation for particle impact to be a significant producer of 
Fe Ka flux. This situation is difficult to reconcile with the observations of hard 
X-ray spectra, or the proposed particle acceleration mechanisms in the accretion 
disk corona. Truncated accretion flows must be externally illuminated by hard 
X-rays in order to produce the Fe Ka line, as proton impact is very inefficient 
in generating line emission. In contrast to the Sun, our conclusion is that, with 
the possible exception for localized regions around magnetic footpoints, particle 
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impact will not be an important contributor to the X-ray emission in accreting 
black holes. 

Subject headings: accretion, accretion disks — line: formation — radiation mech- 
anisms: non-thermal — X-rays: general 

1. Introduction 

An iron Ka emission line at 6.4 keV is commonly observed in the X-ray spectra of active 
galactic nuclei (AGN) and Galactic black hole candidates (GBHCs). The equivalent widths 
(EWs) of the lines are usually ~100-300 eV (Mushotzky, Done & Pounds 1993; Nandra et 
al. 1997; Reeves & Turner 2000), suggesting that a significant column of iron is necessary 
to produce the required emission. This was confirmed with the discovery of the Compton 
refiection hump between 20-30 keV (Pounds et al. 1990; Nandra & Pounds 1994) which 
is a natural outcome of X-ray reprocessing in optically-thick material (Lightman & White 
1988; Guilbert & Rees 1988). In some cases, sensitive observations of the Fe Ka line have 
found them to possess broad low-energy wings, consistent with the profile expected from 
a spinning accretion disk close to the black hole (see the reviews by Fabian et al. (2000) 
and Reynolds & Nowak 2003). These relativistically broadened lines have been seen in both 
AGN (e.g., Tanaka et al. 1995; Nandra et al. 1997, 1999; Fabian et al. 2002) and GBHCs 
(e.g., Martocchia et al. 2002; Miller et al. 2002a,b,c), and provide strong evidence that the 
optically-thick refiector responsible for the shape of the X-ray continuum is the inner region 
of the accretion disk. Thus, there has been considerable attention paid in recent years to 
the Fe Ka line and X-ray refiection features in accreting black holes, as they may probe the 
physics of the fiows only a few Schwarzschild radii from the event horizon. 

The phenomenology outlined above has given rise to a theoretical framework in which 
the X-rays are generated by some process in a corona above the accretion disk (Galeev, 
Rosner & Vaiana 1979; Haardt & Maraschi 1991, 1993; Haardt, Maraschi & GhiseUini 1994). 
Roughly half of the X-ray power illuminates the disk, is reprocessed, and then reemitted as 
a reflection spectrum which carries the Fe Ka line and reflection hump (George & Fabian 
1991; Matt, Perola & Piro 1991). The iron line is produced by fluorescence; that is, incoming 
X-rays (usually in the form of a power-law) with energies greater than 7.1 keV can remove 
a K shell electron from an iron atom or ion. The excited ion can then relax by a radiative 
transition into the K shell, or by the expulsion of an outer shell electron (the non-radiative 
Auger effect) . With this mechanism the fiux in the line is directly proportional to the ionizing 
fiux (Basko 1978; Bai 1979), so that changes in the illuminating continuum should give rise to 
correlated changes in the Fe line fiux. However, long monitoring campaigns of certain bright 
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AGNs have shown that, at least in those sources, the Fe Ka hne exhibits weak variabihty on 
short timescales despite large changes in the continuum (e.g., Chiang et al. 2000; Reynolds 
2000; Vaughan & Edelson 2001). Revisions to the standard disk-corona model are necessary 
to explain this behavior (Reynolds 2001; Fabian & Vaughan 2003; Ballantyne, Vaughan & 
Fabian 2003). 

Fluorescence from photoionization is not the only mechanism able to generate an iron K 
line. Electron bombardment of the accretion disk can also produce line emission from iron if 
the kinetic energy of the particle is greater than 7.1 kcV. This process has been used to explain 
enhanced line emission from the Sun (e.g., Zarro, Dennis & Slater 1992), where electrons 
energized in the solar corona travel down magnetic field lines and strike the photosphere 
(Dennis & Schwartz 1989). As accretion disks hkely also possess such magnetized coronae 
(Merloni & Fabian 2001), it is possible that similar situations would arise and result in Fe Ka 
emission by electron bombardment. 

In accretion disk coronae there exists a bath of thermal electrons with temperature 
kT Ri 100-200 keV (which generates the primary X-ray power-law through Comptoniza- 
tion of the disk UV photons). There may also be a non-negligible amount of non-thermal 
electrons present in the corona. Evidence for these particles has been found in Compton 

Gamma-ray Observatory observations of Cyg X-1 in its soft state (Gierlihski et al. 1999) 
and GRS 1915+105 (Zdziarski et al. 2001). These data, when combined with simultaneous 
hard X-ray observations, show that the X-ray power-law extends to energies > 600 keV. 
Such spectra can be well fit by a model of Comptonization in a hybrid thermal/non-thermal 
plasma (Coppi 1999). At present there are no analogous observations for radio-quiet AGN, 
but it is expected that the disk-corona physics will be similar for the two systems. 

Proton illumination will also produce Fe Ka emission in a manner analogous to electron 
bombardment. This may occur in a geometry where the optically thick accretion disk is 
disrupted into a hot ion flow close to the black hole (e.g., Meyer & Meyer-Hofmeister 1994; 
Rozahska & Czerny 2000a; Manmoto et al. 2000; Manmoto & Kato 2000; Meyer, Liu & 
Meyer-Hofmeister 2000; Spruit & Deufel 2002). This configuration may be relevant to GB- 
HGs in the low/hard state (Esin, McClintock & Narayan 1997; Done & Zycki 1999; Barrio, 
Done & Nayakshin 2003) and low luminosity AGNs (e.g., Quataert et al. 1999). In this sce- 
nario, the accretion energy is envisaged to be dissipated into the protons within the central 
optically thin flow, which then rapidly reach the virial temperature as they can be thermally 
decoupled from the electrons (Shapiro, Lightman & Eardley 1976). The majority of the 
accretion energy is carried into (or advected) into the black hole by the protons with little 
or no observational signatures. Thus, these "advection dominated accretion flows" (ADAFs) 
can be very radiatively inefficient (Ichimaru 1977; Rees et al. 1982; Narayan & Yi 1994). 
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However, if the optically thick disk partially overlaps with the ADAF (for example, in the 
region where evaporation of the disk is taking place), then the hot protons in the ADAF 
will interact with the disk, significantly heating its surface (Spruit & Haardt 2000; Deufel 
& Spruit 2000; Deufel, Dullemond & Spruit 2002) and possibly producing Fe Ka emission. 
Proton impact ionization would also be relevant if a two-temperature corona exists above an 
untruncated accretion disk (Di Matteo et al. 1997). 

If energetic coronal electrons (in an accretion disk-corona model) or virialized protons 
(in a truncated disk- ADAF model) strike the accretion disk they would produce Fe Ka emis- 
sion. If this process is sufficiently efficient it could generate an observable line that is not 
strongly correlated with the continuum (e.g., Petrucci et al. 2002). This paper considers 
Fe Ka production from particle impact onto accretion disks and compares it with the more 
traditional mechanism of photoionization. We consider incident beams consisting of thermal 
or non-thermal electrons and virialized protons. The goal is to evaluate under which circum- 
stances, if any, particle bombardment will be an important contributor to the Fe Ka line. 
The theory needed to calculate the expected line flux is outlined in the next section, and 
Section 3 presents the results. In Section 4 we consider the implications of the computations 
on particle acceleration mechanisms and the two different models of accretion geometry. A 
summary and our conclusions are presented in Section 5. 



2. Mechanisms of Fe Ka Production 

2.1. Electron Bombcirdment 

Consider a beam of electrons incident on the normal of an accretion disk with a spectral 
flux of jVe electrons cm~^ s~^ keV~^ at kinetic energy E. In both AGN and GBHCs the disk 
is hot enough so that both hydrogen and helium are fully ionized. Therefore, each incident 
electron will lose energy by Coulomb interactions with the hydrogen ions as it travels through 
the disk (Emslie 1978): 

dE ^ A, nv , 

where e is the electron charge in e.s.u., and n, v and E are the instantaneous total hydrogen 
number density, electron velocity and energy, respectively. A is the Coulomb logarithm 



where bo is the maximum impact parameter for the scattering event. Equation 1 differs from 
the classic analysis by Spitzer (1962) by ignoring the the correction due to thermal motion 
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of the target electrons (Emslie 1978). This "cold target" approximation holds if £■ » /cT, 
where T is the electron temperature of the accretion disk. Since we are interested in incoming 
electrons with E > x, where x = 7.112 keV is the ionization potential for neutral iron, this 
condition will hold for both AGN {kT - 10 eV) and GBHC {kT - 500 eV) disks. 

As each incident electron decelerates^ from its initial energy Eq it can eject a K shell 
electron from iron as long as it energy remains above x- Thus, the resulting Fe Ka line flux 
in photons cm~^ s~^ is (e.g., Emslie et al. 1986) 



Tei = -(3cu / XEiEo) / aei{E)nF,v{E)dt 

^ Jo Ue=Eo 



. dEo, (3) 

' E=Eq 

where cu is the K fluorescence yield of iron, /3 is the Ka to K/3 branching ratio, aei is the 
K-shell ionization cross-section of iron due to electron impact, and ripe is the number density 
of iron atoms. The inner integral in the above equation counts the number of Fe Ka photons 
produced by an electron with initial energy Eq as it loses energy passing through the disk. 
This then must be integrated over the distribution of initial electron energies to find the 
total Fe Ka line fiux. The factor of one-half in front of the expression accounts for the fact 
that roughly half of the Ka photons produced will be emitted away from the disk surface. 
We assume here, and in the treatment for photoionization described below, that all of the 
line photons emitted toward the surface escape to infinity. 

Equation 1 can be used to rewrite the above expression in a more useful form: 



a^i{E)EdE 



dEo, (4) 



where Apc is the abundance of iron relative to hydrogen in the disk, and we have changed the 
lower limit of the outer integral to the iron K-shell ionization energy (lower energies cannot 
produce Ka lines). Note that this expression is independent of the density distribution of 
the target. 

Of course, electrons with larger energies travel a greater distance before their kinetic 
energy falls below x- Emslie (1978) showed that an electron with initial energy Eq will have 
an energy E = Eq(1 — Qwe'^ANii/ EqY^^ after traversing a hydrogen column A'^h- Therefore, 
the minimum initial energy needed by an electron to just be able to K-shell ionize an iron 
atom after passing through a column Ah is the solution of 

^min - eTre^AATH^^in " = 0. (5) 



^About 10^^ of the energy lost by the electron through scattering is radiated away as bremsstrahlung 
emission in the X-ray band, with the remainder ending up as heat (see, e.g., Tatischeff 2002). This radiation 
is ignored here. 
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Using this relation to replace the infinity in the outer integral of equation 4, will allow us to 
calculate T^i as a function of A^h into the accretion disk. 



2.2. Proton Bombardment 

If protons are the incident particles then the equation for the energy loss rate is modified 
only^ by a factor of mp/mg (Emslie 1978): 

^ = -2vre^^A^, (6) 
dt me E' ^ ^ 

where rUp is the proton mass, nie is the electron mass, and now E is the kinetic energy of the 
impacting proton. Note that protons lose their energy much more rapidly than electrons. As 
before, we have assumed the "cold target" approximation which in this case requires E ^ 
{■mp/me)kT. For AGN disks, this condition implies > 18 keV, while it is > 900 keV 
for the hotter GBHC disks. The virial temperature for a proton a distance r away from a 
black hole is fcT^ir = 156(r/i?s)~^ MeV, where Rs — 2GM/c^ is the Schwarzschild radius 
for a black hole of mass M. Therefore, assuming the protons in the ADAF reach the virial 
temperature, their energies will be greater than 1 MeV if r < 100 i?s and so equation 6 will 
be vahd for both AGN and GBHCs. 

We can now write down the analogous expression to equation 4 for the Fe Ka line flux 
due to proton impact: 



_ 1 (3uA^^ me , TP /p X 



Eo 



api{E)EdE 



dEo, (7) 



where Ve is the spectral proton flux (protons cm~^ s~^ keV~^) normally incident on the disk, 
and (Tpi{E) is the K-shell ionization cross-section of iron due to a proton with kinetic energy 
E. This expression differs notably from equation 4 by a factor of rrie/mp, therefore, unless 
this is compensated by the cross-section, protons will be a much less significant producer of 
Fe Ka photons than electrons. 

After passing through a column A^h, a proton with initial energy Eq will have an energy 
(Emshe 1978) 

1/2 



.^.„(i-Wa(5)^) . ,8) 



^Strictly speaking, A is a function of the reduced mass of the scattering interaction and so wiU change 
with the different incident particle. However, the dependence is logarithmic in nature and so can safely be 
neglected. 
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Therefore, as with the electrons (eq. 5), we can calculate the minimum kinetic energy a 
proton must have to just be able to ionize iron at the bottom of a column A^h^ 

Ernin = ^ + ^T^e^Kr^pl me)Nn. (9) 

Replacing the infinity with £'min(-/VH) in the outer integral in equation 7 gives the total Fe Ka 
fiux produced through a column A^h- 



2.3. Photoionization 

There are many previous studies which have computed the Fe Ka fiuxes expected from 
photoionization (Hatchett & Weaver 1977; Basko 1978, 1979; Bai 1979). In this work, how- 
ever, we will use the simplest expression for the line fiux: 

1 f°° 

•^photo = -NuA^^up J Fe(e)c7photo(e)c?e, (10) 

where F^{e) is the photon flux (ph cm~^ s'^ keV~^) incident on the normal of the accretion 
disk at energy e, and CTphoto is the K-shell photoionization cross-section of iron. As in the 
expression for electron impact production of Fe Ka (eq. 4) , this equation assumes that iron 
is the only metal in the gas. Furthermore, equation 10 is only valid in the optically-thin 
limit. That is, we are ignoring any attenuation of the incident X-rays or the outgoing Fe Ka 
line photons. Therefore, we are limited to columns A^h < 1/(1-2o-t) = 1-25 x 10^^ cm~^, 
where ctt is the Thomson cross-section, and we have assumed a 10% abundance of helium 
(by number). The implications of limiting the calculations to Thomson thin columns are 
considered in the next section. 



2.4. Charge Exchange 

Another particle process that can produce Fe Ka emission is charge exchange collisions. 
These are reactions of the following form: Fe^"*" -|- A ^ Fe*-''^^'^ + A"*", where A is usually 
H, H2, or He. The product Fe ion is often left in an excited state and will produce line 
emission as it relaxes. This collision differs from those described above in that the Fc ion 
is the energetic particle and the hydrogen or helium atom is the target. This process may 
be important for explaining the strong and broad Ka lines from hydrogenic and helium-like 
Fe observed in the Galactic ridge X-ray emission (Tanaka, Miyaji & Hasinger 1999; Tanaka 
2002). In this case, highly ionized Fe ions in the cosmic-ray distribution interacting with 
neutral Galactic gas would produce the needed line emission. 
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For the accretion disk scenario that is considered here, it is unhkely that charge exchange 
will be an important process, as there are very few neutral H and He atoms available for the 
reaction. Other metal ions could act as donors in a collision with Fe, but, for typical cosmic 
abundances, encounters with the correct kinematics would be extremely rare. However, 
charge exchange may possibly be relevant to iron K lines emitted from the molecular torus 
of AGN unification schemes (e.g., Ghisellini, Haardt & Matt 1994; Krolik, Madau & Zycki 
1994) if there is a flux of energetic Fe ions incident on this gas. 



3. Results 

3.1. Electron Impact Versus Photoionization 

To determine how efficient electron impact is in producing an Fe Ka line in comparison 
to photoionization, we calculate the ratio of equations 4 and 10: 



dEo 

■ (11) 



•^photo 27re4AA^H/~-^e(e)(7photo(e)c^e 

Considering this ratio will also help negate the effects of absorption and scattering on the 
Ka line that were ignored in the individual expressions. This value will nevertheless be 
only a upper limit, because in the true Thomson thick case high energy photons can scatter 
down to energies where they can ionize iron and produce a Ka photon. Some fraction of the 
line photons will escape the layer, but to accurately treat this process requires a numerical 
prescription. 

The K-shell photoionization cross-section (cTphoto) was taken from the tabulation by 
Verner & Yakovlev (1995). The empirical fit formula of Hombourger (1998) was used to 
compute the cross-section for K-shcll ionization by electron impact^. These cross-sections 
are plotted between x and 200 keV in Figure 1, where the x-axis is used for both photon 
energy and electron kinetic energy. cXphoto is roughly proportional to e^^ and so falls off much 
more rapidly than however, (Xphoto is significantly greater than cTgi for energies below 
40 keV. 

A power-law photon spectrum was assumed for all the calculations described below: 
F^{e) = Fo(e/x)~'", where F = 1.7, 1.8, 1.9 or 2.0 is the photon index. For each value of F 



•^There is a misprint in the expression for Gr in the paper by Hombourger (1998), so the form quoted by 
Quarles (1976) was used in its place. 
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the spectrum was normalized so that the total energy flux between x and 200 keV at the 
surface of the disk was 1 keV cm~^ s~^. A value of 25 is chosen for the Coulomb logarithm 
A (Brown 1971), but, as will be seen below, the results are insensitive to its value. 

The value of A^h and equation 5 determines the upper-limit to the outer integral in the 
numerator of equation 11. Ka photons produced by electrons with kinetic energies larger 
than this are accounted for in a later bin where the column is large enough for the energy 
of those electrons to pass below x- Since the largest column we consider is 10^"^ cm^^, there 
is a slight error as we do not count any Fe Ka emission produced by electrons that can pass 
through this column. However, equation 5 shows that only electrons with kinetic energies 
greater than ~ 3.5 MeV can still ionize iron beyond a column of 10^"^ cm~^. For the electron 
distributions we consider, there is a neghgible amount of electrons with £■ > 3.5 MeV, so 
this slight error will not impact the results. 

3.1.1. Power-law Afs 

Modeling of the hard X-ray tails in the GBHCs Cyg X-1 and GRS 1915+105 have 
shown that the non-thermal electrons responsible for this emission can be described as a soft 
power-law (Gierhhski et al. 1999; Zdziarski et al. 2001). Therefore, we first consider electron 
distributions of the form jV^ oc [E/x)"'^, where 7 = 1,2,3 or 4 is the number index. As 
with the photon distribution, the incident electron power-laws are normalized such that the 
total energy flux between x and 200 keV is 1 keV cm~^ s~^. In this way, roughly the same 
amount of energy is being deposited into the material by both the electrons and photons. 

Figure 2 shows the results of computing the ratio J-ei/ J'^hoto for the different values of 
r and 7. Each panel shows the results for an individual photon index, while the different 
lines denote the electron number index. In every case considered here the Fe Ka flux from 
electron impact is a small fraction of that produced by photoionization. The maximum ratio 
found is J^ei/ ^photo ~ 0.01 when F = 1.7, 7 = 1 and the Thomson depth (rx = 1.2axA^H) 
is ~ 10~^. The ratio then drops rapidly as the column increases as less and less electrons 
are able to penetrate further into the gas. The only exception is for 7=1. This electron 
spectrum is very flat so that the energy flux per decade actually increases with energy, and so 
there are many high energy electrons that can still ionize Fe Ka at larger column densities. 
The electron K-shell cross-section is close to a constant at high energies (Fig. 1), which helps 
keep the ratio roughly the same as the column increases. However, even with this extreme 
spectrum, electron impact is still an inefficient method to produce Fe Ka photons. 

As 7 is increased, the ratio at all tt drops. This is a result of the shape of the cross- 
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section which peaks at electron energies of ~ 20 keV and not at the the ionization threshold, 
X (Fig- !)• Therefore, when the electron spectrum softens, more electrons have lower energies, 
below the peak in the cross-section, and the Fe Ka flux decreases accordingly. 

There is also a slight change in Jvi/-^photo with the photon- index. The ratio lowers as F 
increases because for softer photon spectra there are more photons at lower energies, which 
is where (Tphoto is larger. This increases the number of Ka photons emitted, and therefore 
decreases the ratio. 

3.1.2. Thermal Ne 

Observations by BeppoSAX have shown that many Seyfert 1 galaxies have high energy 
cutoffs between 100-300 keV (Matt 2001). If the hard X-rays in these AGN are produced 
by thermal Comptonization, then such measurements can constrain the temperature of the 
coronal electrons by kT ^ Ec/2, where is the cutoff energy. We now consider the pro- 
duction of Fe Ka photons by electrons in a thermal distribution, described by a Maxwell- 
Boltzmann distribution with kT = 100, 200, 300 or 400 keV. As before, the electron spectra 
were normalized so that the total energy flux incident on the gas between % and 200 keV 
was 1 keV cm~^ s~^, the same as for the photon spectra. 

Figure 3 plots the Fe Ka flux ratio between electron impact and photoionization for 
electrons in a thermal distribution. Again, the flux of Ka photons produced by the electrons 
is < 1 per cent of that generated by photoionization. An interesting difference between 
the thermal and power-law electron distributions is that the largest flux ratio occurs some 
distance into the layer for a Maxwellian population and not at the surface, as was found with 
the power-law. This is because the number of electrons in the thermal distribution peaks 
at energies ~ kT, and can thus ionize Fe to larger columns. Therefore, the maximum ratio 
moves to greater tt for larger kT. 

Except for the lowest columns, ^ei/-^photo increases with kT. The total number of 
electrons in the distribution grows with kT, and, since the spectrum is normalized so that 
the energy deposited between x and 200 keV is the same for each case, the excess electrons 
are foimd at higher energies. This fact, combined with relative flatness of the K-shell cross- 
section (Fig. 1), explains the small increase in J-^ei/^photo with kT. 
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3.2. Proton Impact Versus Photoionization 

We form the ratio of Fe Ka flux produced by proton impact to that generated by 
photoionization by using equations 7 and 10: 



^photo rup 27re4AiVH /^°^ -F;(e)(Tphoto(e)(ie 

The calculation of the denominator of equation 12 is unchanged from the previous section 
(e.g., A = 25). The K-shell ionization cross-section for proton impact was computed by using 
the empirical formula of Romo- Kroger (1998). This cross-section is plotted along with the 
ones for photoionization and electron bombardment in Figure 4. Unfortunately, the formula 
is only valid for proton kinetic energies less than ~ 50 MeV. Through equation 9, this limits 
the column density that can be probed to be < 3 x 10^^ cm~^. 



3.2.1. Thermal Ve 

This section considers the Fc Ka flux produced by a distribution of protons at the local 
virial temperature, kT^i,-. This situation may be applicable to the case where a relatively 
cold, thin accretion disk overlaps with a hot, ADAF-like flow (Deufel et al. 2002). Due to the 
50 MeV cutoff in the proton impact cross-section, we are limited to proton spectra which do 
not have many particles above this energy. A Maxwell-Boltzmann distribution is assumed 
for the particles, and three cases are considered: kT = 1.56, 3.12 and 6.24 MeV. These values 
correspond to the virial temperature at 100, 50 and 25 -Rs, respectively. The proton spectra 
are normalized so that the total energy flux incident on the gas is 1 keV cm~^ s~^. The 
Fe Ka flux produced by the incident protons will be compared with that from a power-law 
photon spectrum, which is normalized to 1 keV cm~^ s^^ between x ^-nd 200 keV. 

The results are show in Figure 5. For each case, it is found that the protons produce 
a negligible amount of Fe Ka emission as compared to photoionization. This is in part 
due to the normalization conditions. The spectra were normalized so that both the protons 
and the photons deposited the same amount of energy (per unit area, per unit time) into 
the material. As the majority of the protons carried 2-3 orders of magnitude more energy 
than a typical incident photon, there are far fewer protons in the illuminating distribution 
than photons (this should be contrasted with the situation described in § 3.1, where there 
was roughly equal numbers of electrons and photons). However, even if the spectra were 
normalized to the same number flux, the Fe Ka emission produced by protons would still be 
much less than that from photoionization. 
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Of the three values of /cTVir considered, it is the middle one (/cTyir = 3.12 MeV) that 
results in the largest ^pj/J^photo, due to the peak in its energy flux residing closest to the 
maximum cross-section (Fig. 4). On the other hand, it is the kT^i,- = 1.56 MeV spectrum 
that has the greatest ratio at low Tt. This is because this distribution has a large number of 
relatively low energy particles which stop quickly in the gas, and so produce their Ka photons 
close to the surface. The other distributions predominantly have much high energy particles 
which can penetrate through to larger columns. One other effect that is noticeable in Fig. 5 is 
the decrease in J-pi/ J-'photo ^ the photon-index, F, increases. As explained in Sect. 3.1.1, the 
softer photon spectra have more photons at low energies, where the cross-section is larger, 
and will produce more Ka photons, thereby lowering the ratio. 

4. Discussion 

We have shown that particle impact is an inefficient method to produce Fe Ka emission 
as compared to photoionization (see also Basko 1979). In this section, we attempt to place 
this result in the wider context of accretion disk models and particle acceleration mecha- 
nisms. We begin with the electron bombardment results and their relation to models of 
coronae above geometrically thin disks. Proton impact and the two-component disk model 
is discussed at the end of the section. 



4.1. Electrons in Accretion Disk Coronae 

It has been argued for some time that the energy responsible for heating the electrons 
in an accretion disk corona is carried and released by magnetic fields expelled from the 
disk surface (e.g. Galeev et al. 1979; Di Matteo 1998). The magnetic energy (ultimately 
derived from the magnetohydrodynamical (MHD) turbulence in the body of the disk; Miller 
& Stone 2000) would then be liberated by reconnection events with neighboring field lines, 
in a deliberate analogy to flares in the lower solar corona. Some fraction of this energy is 
deposited into electrons which end up in a thermal distribution. The remaining energy would 
then be used to accelerate particles (mostly electrons) to high energies. A key parameter 
is the fraction of energy that is funneled into the non-thermal particles for a given event. 
Unfortunately, as is discussed below, that number is not well constrained. 

X-ray observations have shown that, in the accretion disk environment, the newly ener- 
gized coronal particles give up their energy to the soft X-ray /UV photons from the accretion 
disk. This results in the standard hard X-ray Comptonized power-law that is the hall- 
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mark of these sources. The very hard X-ray tails observed in the spectra of Cyg X-1 and 
GRS 1915+105 show that the non-thermal particles are involved in this process. The results 
of § 3.1 have shown that for electrons to be an important contributor to Fe Ka production 
then ~ 10^-10^ times more energy has to be deposited into the disk by particles then by the 
Comptonized photons. For this scenario to be viable, most of the electrons in the corona 
would have to not lose their energy to Comptonization, implying a very low Compton y pa- 
rameter. However, the observed photon-indices and high energy cutoffs show that |/ ~ 1 in 
most of these systems (Petrucci et al. 2001); therefore, Comptonization is important. Alter- 
natively, due to outflows in the corona, the hard X-rays may be beamed away from the disk 
resulting in a very small radiation flux on the surface (Beloborodov 1999). This scenario is 
also unlikely, as particles would have to be both outflowing (to beam the hard X-rays away 
from the disk) and streaming toward the disk to generate the Fe Ka line. Furthermore, the 
Compton reflection hump seen in many Seyfert Is show that roughly half the observed X-ray 
flux must be intercepted by optically thick material. Perhaps a more fundamental problem 
with a very strong particle flux on the disk surface, is that ^ 99% of the incident energy 
goes into heat and not radiation (Tatischeff 2002). This would lead to a strong expansion of 
the surface layers and perhaps lead to evaporation, as is envisaged in the proton-illuminated 
disks (Deufel et al. 2002), which would certainly affect the observed optically-thick signatures 
such as the Compton-reflection bump. 

On the other hand, the non-thermal particles are hkely to be tied to the magnetic field 
lines in the corona. Thus, there may be a large particle fiux at the position of the magnetic 
footpoints on the surface of the accretion disk. In these very local positions, the electron fiux 
may be much higher than the radiation flux and Fe Ka may be dominated by the electron 
impact process. Indeed, such bright spots are seen on the surface of the Sun as sources 
of hard X-rays (Hudson et al. 1994; Kundu et al. 1995). However, such small regions are 
unlikely to be detected amidst the unresolved accretion disk emission. 

The dominance of the Comptonized hard X-ray power-law in the spectra of accreting 
black holes implies that the acceleration of non-thermal particles in the corona is not highly 
efficient. This conclusion is supported by the above results on Fe Ka production by electron 
impact^. The mechanism by which electrons are accelerated to large kinetic energies in 
accretion disk coronae is not known. Again, one can look to the solar corona for guidance, 
where there are three different processes that are considered to be viable candidates (see 
Miller et al. (1997) and Miller (1998) for recent reviews). The first is stochastic acceleration 



Analogous arguments can be made against conduction from the corona as an important heat source for 
the disk (e.g., Rozaiiska 1999; Rozanska & Czcrny 2000b). In a magnetically dominated corona, the electrons 
are more likely to be Compton cooled before conduction can be effective. 
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by MHD waves in the coronal plasma which are excited by the reconnection events. The 
second is shock acceleration by the plasma jets produced by reconnection (e.g., Blackman 
& Field 1994). However, it is thought that the shocks are unlikely to be strong enough to 
account for most energetic particles (that is, shocks may only "pre-accelerate" ) , although 
they may provide the bulk of the heat that goes into the thermal electrons (Di Matteo 1998). 
The third process is direct acceleration by large-scale electric fields, such as those found in 
reconnecting current sheets. The efficiency of all of these processes is, at this point, largely 
unknown, but numerical reconnection simulations are beginning to show that electric field 
acceleration may be < 10% efficient (Nodes et al. 2003). 

Before leaving non-thermal electrons, we can use some of the basic analytical results of 
electric field acceleration as a consistency check of our conclusions. When thermal electrons 
are subjected to an external electric field they will begin to accelerate along the direction 
of the field, but, for most of the particles, this motion will be limited by Compton drag. 
However, the coUisional drag force on electrons decreases with the velocity, so that particles 
with a large enough speed can "runaway" and be accelerated freely by the electric field. The 
value of the electric field where this critical runaway velocity is equal to the thermal velocity 
Ve of a Maxwellian distribution is called the Dreicer field and is given by (Dreicer 1959, 1960) 



where cUg is the electron plasma frequency. In this situation, the bulk of the electrons in the 
thermal distribution will be accelerated. 

In the context of solar flares, Holman (1985) examined the properties of accelerated 
electrons by sub-Dreicer electric fields^. Since the magnetic field induced by the electron 
beam must be less than the field in the acceleration region (5), Holman (1985) found that 
this limited the acceleration rate (per area) to a maximum of cB/27reL, where L is the length 
of the accelerating current sheet (see also Holman, Kundu & Kane 1989). If the radiation 
flux on the disk is 10^^ erg cm~^ s~^ between 10 eV and 100 keV (Ross, Fabian & Young 
1999; Ballantyne, Ross & Fabian 2001), then, for a F = 1.9 spectrum, this translates to 
~ 10^° erg cm~^ s~^ between x and 200 keV. For a significant Fe Ka flux due to electron 
impact, the energy flux due to electrons must be ~ 10^^ erg cm~^ s"^ (Fig- 2). Converting 
this to a number flux using a 7 = 3 spectrum, and plugging into the acceleration limit 
above gives an estimate of the accelerating magnetic field in the accretion disk coronae, 
B ~ (lO^^L) G. For comparison, the equipartition magnetic field in a disk is ~ 10^ G 



^Super-Dreicer acceleration has been considered by, e.g., Litvinenko (1996) and Craig & Litvinenko (2002), 
where the acceleration of the electrons are limited by transverse magnetic fields. 
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(Merloni & Fabian 2001). Thus, the assumptions of sub-Dreicer acceleration of non-thermal 
particles and significant Fe Ka emission from electron impact implies an extremely strong 
coronal magnetic field. 

4.2. Protons and Truncated Accretion Disks 

Considerable work over the last few years has been done investigating the properties of 
disks which change from being cold and geometrically thin to hot and geometrically thick 
(see references in § 1). One of the greatest difficulties in describing such a model is the 
process by which the truncation of the cold disk occurs. A natural mechanism would be 
heating the disk until its local scale height becomes comparable to the radius, effectively 
evaporating the inner part of the disk. 

Once the two-phase accretion flow is operating, Deufel et al. (2002) showed that protons 
in the hot phase striking the cold disk will be an effective mechanism for heating the disk 
surface. Comptonization of the soft disk photons by this heated layer produced a hard X- 
ray spectrum not too dissimilar from what is observed. In § 3.2 we found that virialized 
protons are a highly inefficient way of producing Fe Ka emission. Thus, external hard X- 
ray radiation will need to be included in the models of truncated accretion flows to explain 
any Fe K emission. In the context of the proton illumination model, it may be possible to 
determine the relative importance of radiation heating to proton heating by examining the 
equivalent width of the Fe Ka line. 

4.3. Other Applications 

The above sections have shown that it is difficult for particle impact to be an efficient 
producer of Fe Ka lines, because most of incoming energy is used to heat the gas and not for 
ionizing iron. Thus, this process will only be important in environments where photons are 
not as efficiently produced as high energy particles. While detailed calculations are beyond 
the scope of this paper, non-thermal radio sources or jet-like flows, where the bulk of the 
energy may be transported by relativistic particles, may be areas where particle impact is 
the dominant source of ionization. The shock front in a gamma-ray burst is another possible 
area of application if there are enough high energy baryons in the relativistic outflow. 
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5. Summciry 

Just like a photon, high energy electrons and protons can eject a K-shell electron from 
iron and produce a Ka line photon if their kinetic energy is > 7.112 keV. In this paper we 
have considered the possible role of particle impact in the generation of Fe Ka emission from 
accreting black holes. Calculations of Fe Ka flux from various columns of gas were performed 
when it was bombarded by a non-thermal power-law of electrons, or a thermal distribution 
of electrons or protons. This flux was then compared with that produced by a power-law of 
radiation under the assumption that all process deposited roughly 1 erg cm~^ s~^ of energy 
into the material. It was found that in every circumstance the Fe Ka flux produced by 
particle impact was < 1% of that produced by photoionization. A more realistic calculation 
which relaxes the optically-thin assumption made here would flnd an even smaller ratio. 

Increasing the electron flux on the accretion disk in order to produce more Ka emission 
was inconsistent with the standard disk-corona model for accretion disk emission. In the 
case of a disk illuminated by virialized protons, as might be found in a truncated accretion 
disk model, we found a very low efficiency for Fe Ka production. Therefore, these type of 
models must consider external X-ray illumination to explain any observed Ka line. 

To conclude, we find that, except for isolated regions around magnetic footpoints, it 
is unlikely that particle impact can significantly contribute to the observed X-ray line (or 
continuum) emission from accreting black holes. Depending on the efficiency of the unknown 
acceleration mechanisms, particles may, however, be a important heating source for accretion 
disks. 

DRB was supported by the Natural Sciences and Engineering Research Council of 
Canada. ACF thanks the Royal Society for support. We thank the anonymous referee 
for useful comments. 



A. Unnormalized Line Fluxes from Pcirticle Impact 

In the interests of generality, we show in Figure 6 the unnormalized Fe Ka line fiuxes (in 
photons cm~^ s^^) obtained from equations 4 and 7 for the power-law and thermal spectral 
shapes used in § 3. The standard values of a; = 0.342 and P — 0.822 from Bambynek et al. 
(1972) were used, and the solar Fe abundance of Ap^ — 2.82 x 10~^ from Holweger (2001) 
was assumed. As the Fe abundance is just a multiplicative factor, the results can be easily 
scaled to other metallicities. Finally, the particle spectra were normalized as in § 3; that is, 
the electron spectra has a energy flux between 7.112 keV and 200 keV of 1 keV cm~^ s~^. 
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while the protons have a total energy flux of 1 keV cm ^ s ^. Again, the results in Fig. 6 
can be easily scaled to other impacting fluxes. 
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Fig. 1. — K-shell ionization cross-section of iron as a function of energy for two different 
ionization processes. The solid line shows the cross-section for photoionization using the fits 
of Verner & Yakovlev (1995). The dotted line denotes the empirical expression for electron 
impact cross-section found by Hombourger (1998) (in this case the x-axis plots electron 
kinetic energy). 
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Fig. 2. — Ratio of Fe Ka flux from electron impact J^ei to Fe Ka flux from photoionization 
•^photo- Both incident spectra were power-laws with the radiation having a photon index F 
and the electrons a number index 7. Both spectra were normalized so that the energy flux 
between x and 200 keV was 1 keV cm~^ s~^. In each panel J-'ei/J^photo is plotted for 7 = 1 
(solid line), 2 (dotted line), 3 (short-dashed line), and 4 (long-dashed line). In all cases the 
Fe Ka line flux produced by electron impact is negligible compared to the one produced by 
photoionization. 



-24- 




0.001 0.01 0.1 0.001 0.01 0.1 



Thomson Depth Thomson Depth 

Fig. 3. — Ratio of Fe Ka flux from electron impact J^ei to Fe Ka flux from photoionization 
•^photo- The incident radiation spectrum was a power- law with a photon index F. The 
incident electron spectrum was a Maxwell- Boltzmann distribution of temperature kT. Both 
spectra were normalized so that the energy flux between x and 200 kcV was 1 keV cm~^ s~^. 
In each panel J-'ei/^photo is plotted for kT = 100 keV (solid line), 200 keV (dotted line), 
300 keV (short-dashed line), and 400 keV (long-dashed line). In all cases the Fe Ka line flux 
produced by electron impact is negligible compared to the one produced by photoionization. 
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Fig. 4. — Comparison of the Fe K-shell ionization cross-section due to proton impact (dot- 
dashed hne) with those for electron impact (dotted hne) and photoionization (sold line). The 
proton impact cross-section was taken from the empirical fit of Romo-Kroger (1998). This 
fit is only valid to proton kinetic energies of ~50 MeV, so the cross-section is truncated at 
that point. 
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Fig. 5. — Ratio of Fc Ka flux from proton impact J-pi to Fe Ka flux from photoionization 
•^photo- The incident radiation spectrum was a power-law with a photon index F. The incident 
proton spectrum was a Maxwell-Boltzmann distribution of temperature kT. The radiation 
spectrum was normahzed so that the energy flux between x and 200 keV was 1 keV cm"^ s"^. 
The proton spectrum was normahzed so that the total energy flux was 1 keV cm~^ s~^. In 
each panel J^pi /J^photo is plotted for kT — 1.56 MeV (solid line) corresponding to the virial 
temperature at 100 Rs, 3.12 MeV (dotted hne; 50 Rs), and 6.24 MeV (short-dashed line; 
25 Rs). In all cases the Fe Ka hne flux produced by proton impact is neghgible compared to 
the one produced by photoionization. The cutoff in the proton ionization cross-section (see 
Fig. 4) limits the calculations to low Thomson depth and low virial temperatures. 
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Fig. 6. — Iron Ka flux in photons cm"^ produced by particle impact as a function of 
hydrogen column density. A solar abundance of iron has been assumed, (a) The line flux 
produced by electron impact when the particles have a power-law spectrum. The energy 
flux between 7.112 keV and 200 keV was 1 keV cm~^ s"^. The line styles are the same as 
in Fig. 2. (b) The hne flux produced by electron impact when the particles have a thermal 
spectrum. The energy flux between 7.112 keV and 200 keV was 1 keV cm""^ s--*^. The line 
styles are the same as in Fig. 3. (c) The line flux produced by proton impact when the 
particles have a thermal spectrum. The total energy flux was 1 kcV cm-^ s-^. The line 
styles are the same as in Fig. 5, and the 50 McV cutoff in the proton ionization cross-section 
limited the calculation to lower column densities. 



